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[i]  Electron  acceleration  inside  the  Earth’s  magnetosphere  is  required  to  explain  increases 
in  the  ~MeV  radiation  belt  electron  flux  during  magnetically  disturbed  periods.  Recent 
studies  show  that  electron  acceleration  by  whistler  mode  chorus  waves  becomes  most 
efficient  just  outside  the  plasmapause,  near  L  =  4.5,  where  peaks  in  the  electron  phase 
space  density  are  observed.  We  present  CRRES  data  on  the  spatial  distribution  of  chorus 
emissions  during  active  conditions.  The  wave  data  are  used  to  calculate  the  pitch  angle 
and  energy  diffusion  rates  in  three  magnetic  local  time  (MLT)  sectors  and  to  obtain  a 
timescale  for  acceleration.  We  show  that  chorus  emissions  in  the  prenoon  sector  accelerate 
electrons  most  efficiently  at  latitudes  above  15°  for  equatorial  pitch  angles  between  20° 
and  60°.  As  electrons  drift  around  the  Earth,  they  are  scattered  to  large  pitch  angles  and 
further  accelerated  by  chorus  on  the  nightside  in  the  equatorial  region.  The  timescale 
to  accelerate  electrons  by  whistler  mode  chorus  and  increase  the  flux  at  1  MeV  by  an  order 
of  magnitude  is  approximately  1  day,  in  agreement  with  satellite  observations  during  the 
recovery  phase  of  storms.  During  wave  acceleration  the  electrons  undergo  many  drift 
orbits  and  the  resulting  pitch  angle  distributions  are  energy-dependent.  Chorus  scattering 
should  produce  pitch  angle  distributions  that  are  either  flat-topped  or  butterfly-shaped.  The 
results  provide  strong  support  for  the  wave  acceleration  theory. 
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1.  Introduction 

[2]  The  flux  of  energetic  electrons  in  the  Earth’s  outer 
radiation  belt  can  vary  by  1000  fold  or  more  during 
magnetic  storms  and  other  disturbances  at  the  Earth  driven 
by  the  Sun  [Baker  et  al.,  1986,  1997,  1998a,  1998b].  These 
variations  are  due  to  a  competition  between  acceleration  and 
loss  processes  acting  within  the  Earth’s  magnetic  field  [Li  et 
al. ,  1997;  Baker  et  al. ,  1989].  The  flux  variations  are 
important  since  they  have  been  associated  with  damage 
and  loss  of  satellites  at  geostationary  orbit  [e.g.,  Baker  et  al.y 
1998c;  Baker ,  2001;  Horne ,  2003].  Therefore  to  develop 
physics-based  models  to  analyze  and  predict  the  Earth’s 
radiation  environment,  it  is  important  to  understand  the 
acceleration  and  loss  processes. 
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[3]  Several  theories  have  been  suggested  to  account  for 
the  acceleration  (see  the  reviews  by  Li  and  Temerin  [2001], 
Friedel  et  al.  [2002],  and  Horne  [2002]).  For  very  rapid  flux 
increases  on  timescales  of  minutes,  shock  acceleration  has 
been  suggested,  since  it  can  account  for  the  very  fast 
timescales  [Hudson  el  al.y  1997],  However,  these  events 
are  relatively  unusual.  Far  more  common  arc  flux  increases 
and  decreases  associated  with  magnetic  storms  [Reeves  et 
al.t  2003]  and  periods  of  prolonged  substorms  [Meredith  et 
al.,  2002a,  2003a],  which  have  timescales  for  flux  increase 
of  the  order  of  many  hours  to  days. 

[4]  Betatron  acceleration  as  a  result  of  inward  radial 
diffusion  and  conservation  of  the  first  adiabatic  invariant 
is  an  important  mechanism  that  operates  inside  the  radiation 
belts  [Schulz  and  Lanzerottiy  1974].  This  mechanism  can  be 
enhanced  by  UTF  waves  [Elkington  et  al.,  1999].  However, 
there  is  increasing  evidence  from  peaks  in  the  phase  space 
density  [Brauligam  and  Albert,  2000;  Green  and  Kivelson, 
2004]  and  pitch  angle  distributions  [Horne  et  al.,  2003a] 
that  local  acceleration  by  another  mechanism  is  required  in 
the  heart  of  the  outer  radiation  belt  near  L  ~  4.5.  One  of  the 
most  promising  candidates  is  local  acceleration  by  whistler 
mode  waves  through  Doppler-shifted  cyclotron  resonance 
[Horne  and  Thorne ,  1998;  Summers  et  al.,  1998].  The 
model  proposed  is  that  substorms  and  inward  radial  diffu¬ 
sion  provide  a  seed  population  of  low-energy  electrons.  As 
they  are  injected  to  higher  magnetic  field  strength,  the 
distribution  develops  an  anisotropy  through  conservation 
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of  the  first  two  adiabatic  invariants  and  excites  whistler 
mode  waves.  These  waves  accelerate  the  particles  by 
transferring  energy  from  low-energy  electrons  at  small  pitch 
angles  to  high-energy  electrons  at  large  pitch  angles  through 
Doppler-shifted  cyclotron  resonance.  Indeed,  several  studies 
have  shown  that  relativistic  electron  flux  enhancements  are 
associated  with  prolonged  substorm  activity,  an  elevated 
flux  of  seed  electrons,  and  prolonged  periods  of  enhanced 
chorus  amplitudes  lasting  for  the  order  of  several  days 
[Meredith  et  al. ,  2002a,  2003a;  Miyoshi  et  a!.,  2003].  A 
recent  statistical  analysis  of  magnetically  disturbed  periods 
suggests  that  radial  diffusion  is  efficient  in  the  outer 
magnetosphere,  whereas  local  acceleration  may  dominate 
in  the  inner  magnetosphere  [O’Brien  et  al ,  2003]. 

[5]  In  order  for  whistler  mode  acceleration  to  be  viable, 
waves  must  be  able  to  accelerate  the  electrons  on  timescales 
consistent  with  observations.  While  previous  estimates  have 
generally  been  favorable,  they  have  assumed  field-aligned 
waves  [Summers  et  al,  2002,  2004]  and  a  uniform  distri¬ 
bution  of  wave  power  in  magnetic  local  time  (MLT)  [Horne 
et  al,  2003b].  Another  factor  that  is  critical  to  the  efficiency 
of  wave  acceleration  is  the  ratio  of  the  electron  plasma 
frequency  to  the  electron  gyro  frequency  fpjfce  [Summers  et 
al,  1998;  Home  et  al,  2003b].  Whistler  mode  acceleration 
is  more  efficient  in  regions  of  low  density,  since  this 
increases  the  phase  velocity  of  the  waves  for  the  dominant 
cyclotron  resonance.  Observations  show  that  whistler  mode 
wave  power  and  fpJfce  vary  considerably  with  magnetic 
activity,  L,  MLT,  and  magnetic  latitude  [Meredith  et  al, 
2003c].  The  purpose  of  this  paper  is  to  take  these  variations 
into  account  by  analyzing  data  from  the  CRJRES  satellite  in 
different  MLT  sectors  and  constructing  models  to  calculate 
the  timescale  for  electron  acceleration  on  a  global  scale.  To 
accomplish  this,  we  have  developed  new'  pitch  angle  and 
energy  diffusion  codes  that  are  not  restricted  to  the  high 
plasma  density  approximation  made  in  previous  codes  but 
can  also  treat  diffusion  in  the  low-density  limit  [Glauert  and 
Horne,  2005;  Albert ,  2005].  We  find  that  wave  acceleration 
is  very  dependent  on  magnetic  local  time  and  latitude,  but 
the  resulting  timescale  should  be  sufficient  to  increase  the 
~MeV  electron  flux  by  an  order  of  magnitude  at  L  «  4.5 
within  24  hours,  consistent  with  observations. 

2.  Data  Distribution  in  MLT 

[6]  Using  data  from  the  Combined  Release  and  Radiation 
Effects  Satellite  (CRRES),  we  constructed  a  database  to 
study  the  spatial  distribution  of  whistler  mode  wave  power 
and  the  ratio  fpJfce .  CRRES  operated  in  an  elliptical 
geosynchronous  transfer  orbit  with  a  perigee  of  305  km, 
apogee  of  35,768  km,  and  inclination  of  18°.  The  orbital 
period  was  approximately  10  hours,  and  the  mission  lasted 
for  almost  15  months.  Thus  during  the  mission  lifetime,  the 
satellite  sampled  all  MLT  through  the  radiation  belts.  The 
data  used  in  this  study  were  taken  from  the  Plasma  Wave 
Experiment,  which  measured  electric  fields  between  5.6  Hz 
and  400  kHz  [Anderson  et  al,  1992],  and  the  fluxgate 
magnetometer  [Singer  et  al,  1992].  Since  pitch  angle  and 
energy  diffusion  rates  scale  as  the  power  of  the  wave 
magnetic  field,  the  electric  field  spectral  intensity,  Se ,  was 
converted  into  magnetic  field  spectral  intensity,  SB,  by  using 
Faraday’s  law  and  the  cold  plasma  dispersion  relation  for 


wave  propagation  parallel  to  the  ambient  magnetic  field  B0 
[Stix,  1 992,  p.  26],  given  by 

where  c  is  the  speed  of  light  and  /  is  the  wave  frequency. 
The  plasma  frequency  was  determined  from  observations  of 
electrostatic  weaves  at  the  local  upper  hybrid  frequency  funR 
and  from  the  lower- frequency  cutoff  of  electromagnetic 
continuum  radiation,  as  described  by  Meredith  et  al 
[2002a].  The  magnetic  field  wave  amplitude  was  obtained 
by  integrating  the  average  wave  spectral  density  (pT2  Hz"  *) 
over  frequency  and  taking  the  square  root.  The  data  were 
sampled  and  binned  in  steps  of  0.1Z,.  The  L  shell  and 
magnetic  latitude  \m  were  obtained  from  the  Olsen- Pfitzer 
magnetic  field  model  [Olsen  and  Pfitzer,  1977]. 

[7]  The  distribution  of  chorus  wave  amplitudes  with  jptJ 
fce  have  been  presented  for  different  levels  of  substorm 
activity  before  [Meredith  et  al,  2003c].  The  previous 
studies  used  the  level  of  substorm  activity,  as  measured 
by  the  AE  index,  and  selected  data  for  AE  >  300  nT  as  the 
highest  level  of  substorm  activity.  They  also  limited  the 
frequency  band  to  0.\fce  </<  0.5  fce.  However,  since 
diffusion  rates  are  sensitive  to  the  peak  wave  power  and, 
as  we  shall  show,  are  sensitive  to  lower  normalized  fre¬ 
quencies  at  high  latitudes,  we  have  reanalyzed  the  CRRES 
data  for  AE  >  500  nT,  andfu/R  </<  0.5  fce ,  where  fum  is 
the  lower  hybrid  resonance  frequency.  Therefore  in  this 
paper  we  define  lower-band  chorus  as  chorus  in  the  fre- 
quency  range  fum  </<  0.5  fce  and  not  0.1  fce  </<  0.5  /„ 
used  in  previous  studies. 

[s]  The  average  lower-band  chorus  amplitudes  during 
active  conditions  are  shown  as  a  function  of  L  and  MLT 
for  the  equatorial  and  midlatitude  regions  in  the  top  panels 
of  Figure  L  The  analysis  is  restricted  to  lower-band  choais 
since  these  waves  resonate  with  ~  MeV  electrons  and  since 
lower-band  chorus  is  generally  stronger  than  upper-band 
chorus  [Meredith  et  al,  2001,  2002b].  In  the  equatorial 
region  (top  left),  defined  here  as  waves  between  --15°  < 
\m  <  15°,  chorus  wave  amplitudes  are  enhanced  from  just 
before  midnight,  through  dawn,  to  about  1200  MLT.  Wave 
amplitudes  are  enhanced  for  3  <  L  <  6  and  can  exceed 
100  pT.  Unfortunately,  there  are  very  few  samples  in  the 
0900-1200  MLT  sector.  The  region  of  enhanced  wave 
amplitudes  also  corresponds  to  regions  where  fptJfce  is 
typically  less  than  4  (bottom  left),  indicating  that  the  region 
from  midnight  through  dawn  should  contribute  to  wave 
acceleration.  In  comparison  to  equatorial  chorus,  midlati¬ 
tude  chorus  (between  15°  <  |\m|  <  30°)  is  enhanced  on  the 
dayside  between  0600  and  1500  MLT.  The  emissions  before 
1200  MLT  are  associated  with  low  values  of fpJfce,  typically 
less  than  3,  and  thus  should  also  contribute  to  wave 
acceleration.  However,  at  night  there  is  very  little  wave 
power  at  high  latitudes  even  though  the  ratio  fpJfce  is 
smaller.  The  data  show  that  there  are  substantial  differences 
in  wave  amplitudes  and  fpJfce  with  MLT,  and  thus  there 
should  be  substantial  differences  in  the  pitch  angle  and 
energy  diffusion  rates. 

[9]  To  take  into  account  MLT  effects  on  the  diffusion 
rates,  we  have  averaged  the  CRRES  data  into  three  local 
time  sectors,  where  in  each  sector  the  data  exhibit  similar 
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Figure  1.  Magnetic  local  time  (MLT)  dependence  of  (top)  chorus  wave  amplitudes  and  (bottom)  the 
electron  plasma  to  gyro  frequencies  fpJfce-  The  data  are  for  high  substorm  activity  AE  >  500  nT,  for  the 
(left)  equatorial  and  (right)  midlatitude  regions.  The  insets  show  the  time  spent  sampling  each  0.1L  bin. 


variations  in  wave  power  and  latitude.  The  MLT  sectors  are 
referred  to  as  night  (2300-0600  MLT),  prenoon  (0600- 
1200  MLT),  and  afternoon  (1200-1500  MLT).  The  remain¬ 
ing  MLT  sector  is  not  included,  since  the  wave  power  is  low 
and  fpJfCe  is  high  so  that  no  significant  acceleration  is 
expected.  The  data  are  restricted  to  L  =  4.5  ±0.5L, 
corresponding  to  the  approximate  peak  in  electron  phase 
space  density  during  the  recovery  period  of  magnetic  storms 
[Green  and  Kivelson ,  2004].  All  the  data  are  for  a  high  level 
of  substorm  activity,  AE  >  500  nT,  and  data  for  southern 
latitudes  has  been  averaged  with  northern  latitudes  to 
improve  the  statistics.  Figure  2  shows  the  results  for  the 
median  (solid),  upper,  and  lower  (dashed)  quartiles  of  the 
average  wave  amplitude  for  fum  </<  0.5  fce  (top)  and  ratio 
fpjfce  (bottom). 

[io]  Chorus  wave  amplitudes  peak  around  50  pT  for  \m  < 
15°  in  the  night  MLT  sector  and  decrease  with  increasing 
latitude  (top  left).  On  the  other  hand,  wave  amplitudes 
increase  with  latitude  in  the  prenoon  sector,  reaching  a 
maximum  around  \m  «  25°.  The  prenoon  sector  has  the 
largest  wave  amplitudes  of  all  the  MLT  shown.  In  the 
afternoon  sector,  wave  amplitudes  also  tend  to  increase 
with  increasing  latitude  but  are  generally  smaller  than 
prenoon.  The  ratio  fptJfce  generally  decreases  with  increas¬ 
ing  latitude  in  each  MLT  sector.  However,  fpJfce  is  smaller 
in  the  night  sector.  The  data  for  fpJfcc  are  more  variable  in 
the  prenoon  and  afternoon  sectors,  probably  due  to  a  time 


delay  between  the  increase  in  AE  above  500  nT  and  the 
motion  of  the  high-density  plasmapause. 

[it]  Using  the  data,  we  have  constructed  three  models  for 
the  variation  in  wave  amplitude  and  plasma  density 
corresponding  to  each  MLT  sector  at  L  =  4.5.  We  have 
restricted  the  magnetic  field  variations  to  a  dipole  and  kept 
the  plasma  density  constant  with  latitude.  The  wave  power 
in  each  model  is  kept  constant  to  reduce  the  complexity  of 
the  calculations  and  to  help  physical  interpretation.  The 
models  are  overlaid  in  Figure  2  (long  dashed  lines). 

[12]  For  the  night  model,  fpJfce  ~  3.43  at  the  magnetic 
equator  and  the  average  wave  magnetic  field  amplitude  is 
Bw  =  50  pT  over  a  latitude  range  |\m|  <  15°.  The  wave 
amplitude  is  zero  outside  this  range.  For  the  prenoon  model, 
fpjfce  =  4.0  at  the  equator  and  Bw  =  100  pT  for  15°  <  \\m\  < 
35°.  This  gives  a  reasonable  variation  over  the  latitude 
region  where  the  wave  power  peaks,  for  example,  fpJfce  = 
1.8  at  \m  =  25°.  For  the  afternoon  mod  el,  fptJfce  =  6.72  at  the 
equator  and  Bw  =  50  pT  for  10°  <  |\m|  <  35°. 

3.  Diffusion  Model 

[13]  Following  previous  methods  [Lyons  et  al ,  1972; 
Lyons ,  1974;  Albert ,  1999],  we  calculate  the  bounce- 
averaged  diffusion  coefficients  for  equatorial  pitch  angle 

mixed  pitch  angle  momentum  {D(x  p),  and  mo¬ 
mentum  (Dpp)  in  the  fully  relativistic  limit.  The  bounce- 
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Chorus  Band:  (  f  (  0.5f 

lhr  N  '  ce 

L  Coverage:  4.5  ±  0.5;  Magnetic  Activity:  AE  )  500  nT 

23:00  -  06:00  MLT  06:00  -  12:00  MLT  12:00  ->  15:00  MLT 


Magnetic  latitude 


Magnetic  latitude 


23:00  -  06:00  MLT  06:00  -  12:00  MLT 


12:00  -  15:00  MLT 


-  median 

.  upper  and  lower  quartiles 


model  values 


Figure  2.  (top)  Average  chorus  wave  amplitudes  and  (bottom)  ratio  fpJfce  for  different  magnetic  local 
time  sectors,  (left)  night,  (middle)  prenoon,  and  (right)  afternoon,  for  AE  >  500  nT  at  L  =  4.5.  Data  for 
southern  and  northern  latitudes  are  averaged  together.  The  median,  upper,  and  lower  quartiles  are  shown 
by  the  solid  and  short  dashed  lines,  respectively.  The  model  values  are  shown  by  the  long  dashed  lines. 


averaged  diffusion  coefficients  are  obtained  by  bounce- 
averaging  the  relativistic  quasi-linear  diffusion  equation, 
given  by  [Lyons  and  Williams ,  1984,  p.  140] 


dfo 

dt 


V.(D.V/0) 

_J _ d_ 

p 2  sin  ot<9a 
J  d  2( 
+  P2dPP  { 


/?  sin 


a  -dJ±  +  D  dA 

poo.  op 


D 

paPdK  pp  dp 


(2) 


where  a  is  the  particle  pitch  angle  (tan  a  -  pJp\\),fo(p,  0  is 
the  zero-order  spatially  uniform  particle  distribution  func¬ 
tion  for  the  species  cr,  with  charge  qa,  and  rest  mass  m0a, 
and  Da0o  and  Dpp  are  the  local  diffusion  coefficients 
defined  as 


„  />2<(Aa)2> 

Uaa  ~  2A( 

(3) 

„  f(AotAp) 

v  _  2A( 

(4) 

n  <(Ap)2} 

2A 1  ■ 

(5) 

It  is  important  to  note  that  the  diffusion  rates  have 
dimensions  of  momentum2  s-1.  However,  in  the  results 
below,  we  divide  by  p 2  to  give  the  diffusion  rates  in  s~\ 

The  energy  ( DEE )  and  mixed  pitch  angle  energy  (D^) 


diffusion  rates  are  obtained  from  the  momentum  diffusion 
rates  and  divided  by  E2  to  obtain  the  same  units  of  s-1. 

[m]  While  previous  work  uses  the  high-density  approx¬ 
imation  for  whistler  mode  dispersion,  we  use  the  full 
dispersion  relation.  This  enables  us  to  address  diffusion  in 
low  and  high  plasma  densities  but  increases  the  mathemat¬ 
ical  complexity  significantly.  The  results  have  been 
obtained  and  checked  by  two  new  and  independent  codes. 
Both  use  the  full  cold  plasma  dispersion  relation,  but  one 
(PADIE)  solves  a  high-order  polynomial  to  obtain  all 
possible  resonant  frequencies  and  then  selects  those 
corresponding  to  the  whistler  branch  [ Glauerl  and  Horne , 
2005],  while  the  other  uses  the  properties  of  the  whistler 
mode  to  find  the  resonances  within  the  specified  range  of 
frequency  and  wavenormal  angles  [Albert,  2005]. 

[is]  To  calculate  the  diffusion  coefficients,  the  total  wave 
power,  frequency  distribution,  and  angular  distribution  of 
wave  power  must  be  specified.  The  frequency  distribution  is 
assumed  to  be  Gaussian,  given  by 


£2(oj)  =  A2  exp 


'W  —  UJm 
<  6u) 


) 


2 


(6) 


where  £2(w)  is  the  power  spectral  density  of  the  wave 
magnetic  field  (in  nT2  Hz-1),  and  8 to  are  the  frequency  of 
maximum  wave  power  and  bandwidth,  respectively  (in 
radians  s-1),  and 


I  Bw 


6lj  y/tt 


<7> 
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Bw  is  in  units  of  nT.  The  wave  spectrum  is  also  bounded  by 
upper  ujw£.  and  lower  u>/c  frequency  limits  and  is  zero  outside 
these  limits.  We  also  assume  that  the  wave  power  is  spread 
over  a  Gaussian  distribution  of  wave  normal  angles  ^  given 
by 


gUx)  =  exp 


’x-xmy 


(8) 


where  X-  tan  Xw  is  the  angular  width  and  Xm  is  the  peak. 
The  derivation  of  the  diffusion  coefficients,  and  methods  of 
computation,  are  given  in  full  elsewhere  [Glauert  and 
Horne ,  2005;  Albert ,  2005]. 

[i6]  Since  growth  rates  for  whistler  mode  waves  maxi¬ 
mize  for  parallel  propagation,  for  each  model  we  assume 
that  the  wave  power  peaks  in  the  field-aligned  direction, 
Xm  =  0,  with  an  angular  distribution  of  =  30°  ( Xw  - 
0.577).  At  each  frequency  the  wave  power  is  integrated 
from  X  =  0  to  X  ~  1  (or  up  to  the  resonance  cone  angle  if 
that  is  smaller).  Since  wave  propagation  at  an  angle  to  the 
magnetic  field  is  included,  we  calculate  the  diffusion  rates 
for  Landau  (/?  =  0)  and  ±5  cyclotron  harmonic  resonances. 

4.  Bounce-Averaged  Diffusion  Rates 

[n]  The  bounce-averaged  pitch  angle  and  energy  diffu¬ 
sion  rates  for  the  night  model  are  shown  in  Figure  3.  For  the 
night  model  we  assume  that  the  wave  amplitude  peaks  at 
wm  -  0.35£2l,  and  has  a  bandwidth  6oj  =  0.1 5Qey  typical  of 
lower-band  chorus  emissions  [e.g.,  Meredith  et  al.,  2001]. 
The  lower  and  upper  cutoff  frequencies  are  set  to  1 0fc  = 
0.0512e  and  ur/c  =  0.65f2e,  respectively.  The  momentum 
diffusion  rates  (Dpp)  have  been  converted  into  energy 
diffusion  rates  (Dee)' 

[is]  The  very  narrow  peaks  in  the  diffusion  rates  near  80° 
correspond  to  diffusion  by  individual  resonances,  particu¬ 
larly  the  n-  0  Landau  and  n  =  ±  1  cyclotron  resonances.  The 
higher-order  resonances  occur  at  smaller  pitch  angles  and 
add  in  such  a  way  as  to  give  a  smoother  profile.  The  largest 
pitch  angle  diffusion  rates  occur  at  energies  of  ~30  keV.  At 
low  energies,  pitch  angle  diffusion  maximizes  at  small  pitch 
angles  and  thus  electron  scattering  into  the  loss  cone  is  very 
effective.  Energy  diffusion  maximizes  for  <  60°  (bot¬ 
tom)  and  thus  for  any  moderate  loss  cone  anisotropy,  30  keV 
electrons  contribute  the  most  to  whistler  mode  wave  growth 
as  they  are  scattered  to  smaller  pitch  angles.  At  higher 
energies,  pitch  angle  diffusion  is  reduced  at  small  pitch 
angles  and  peaks  between  60°  and  80°.  This  indicates  that 
pitch  angle  diffusion  is  most  effective  at  large  angles  but 
that  the  particles  remain  trapped  in  the  magnetosphere, 
typically  for  periods  of  several  days  for  E>  300  keV.  These 
trapped  particles  can  undergo  significant  energy  diffusion. 
For  example,  assuming  a  separation  of  variables,  pitch  angle 
and  energy,  in  the  diffusion  equation,  the  timescale  t e  for 
electrons  to  diffuse  in  energy  by  a  factor  of  e  (2.718)  by 
pure  energy  diffusion  can  be  estimated  from  the  inverse  of 
the  energy  diffusion  coefficient.  Thus  the  energy  diffusion 
timescale  for  300  keV  electrons  is  less  than  13  hours. 

[19]  The  net  particle  diffusion  flux  is  toward  lower  phase 
space  density.  At  large  pitch  angles  and  for  a  particle 
distribution  function  that  has  a  modest  anisotropy,  the 


Night  model 


Pitch  Angle  aeq° 


Figure  3.  Bounce-averaged  pitch  angle  (D^O  and  energy 
(DEE)  diffusion  rates  for  whistler  mode  chorus  waves  as  a 
function  of  equatorial  pitch  angle  at^,  for  energies  of  30, 
100,  300,  and  1000  keV,  for  the  night  model. 


resonant  diffusion  surfaces  along  which  the  particles  diffuse 
cut  across  the  contours  of  the  distribution  function  [e.g., 
Horne  and  Thorne ,  2003,  Figure  3].  Under  these  conditions, 
pitch  angle  diffusion  at  large  pitch  angles  results  in  a  net 
flux  of  particles  toward  higher  energies  and  thus  energy 
gain.  Thus  in  Figure  3,  there  should  be  a  net  acceleration 
of  particles  at  pitch  angles  more  than  ~60°.  In  the  absence 
of  other  source  and  loss  processes,  the  diffusion  should 
continue  until  the  contours  of  the  particle  distribution 
function  match  the  resonant  diffusion  surfaces.  The  net 
result  is  a  transfer  of  energy  from  the  large  number  of 
low-energy  particles  at  ~30  keV,  into  the  waves,  and 
acceleration  of  a  smaller  number  of  >100  keV  electrons  at 
large  pitch  angles. 

[20]  In  comparison  to  the  night  model,  the  diffusion  rates 
for  the  prenoon  model  (Figure  4)  show  a  different  behavior. 
In  the  prenoon  model  the  frequency  of  maximum  wave 
power  at  the  equator  is  larger,  =  0.55f2e,  in  order  to 
obtain  the  appropriate  frequency  band  for  chorus  at  higher 
latitudes  |\„,|  >  15°  where  wave  amplitudes  maximize.  The 
other  parameters  are  =  0.1 50e,  ufc  =  0.05f2e,  and  u)Mf  = 
0.95Qe.  In  this  model  all  the  diffusion  rates  drop  rapidly  for 
pitch  angles  greater  than  about  55°.  This  is  due  to  our  model 
assumption  that  there  is  no  wave  power  at  latitudes  less  than 
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Pre-noon  model 


Figure  4.  Bounce-averaged  pitch  angle  (Daa)  and  energy 
(Dee)  diffusion  rates  for  whistler  mode  chorus  waves  as  a 
function  of  equatorial  pitch  angle  aeq ,  for  energies  of  30, 
100,  300,  and  1000  keV,  for  the  prenoon  model. 

15°,  and  thus  when  the  diffusion  rates  are  mapped  back  to 
the  equator  there  is  no  scattering  of  the  equatorial  distribu¬ 
tion  function  at  pitch  angles  greater  than  60°.  Significant 
pitch  angle  scattering  into  the  loss  cone  is  possible  for 
energies  from  30  to  100  keV,  with  a  timescale  comparable  to 
hour,  but  the  losses  become  small  for  energies  above 
300  keV.  In  particular,  at  1  MeV  loss  is  insignificant, 
suggesting  that  once  particles  are  accelerated  to  this  energy 
they  remain  trapped  for  ~10  days.  The  largest  energy 
diffusion  (Figure  4,  bottom)  occurs  for  ~30  keV.  The 
estimated  timescale  for  acceleration  300  keV  electrons  is 
less  than  5  hours  but  decreases  with  increasing  energy  to 
about  1  day  at  1  MeV.  The  acceleration  is  only  important  for 
pitch  angles  between  about  30°  and  60°  at  these  higher 
energies.  Therefore  the  prenoon  model  alone  does  not 
produce  acceleration  all  the  way  to  90°.  In  reality,  there  is 
wave  power  at  low  latitudes  in  the  prenoon  sector  and  so 
there  will  be  some  energy  diffusion  at  large  pitch  angles.  In 
addition,  pitch  angle  diffusion  is  usually  much  larger  than 
energy  diffusion,  and  so,  even  though  the  power  may  be 
small  near  the  equator,  pitch  angle  diffusion  will  tend  to 
smooth  out  any  accelerated  electrons  so  they  should  be  seen 
near  the  equator. 

[21]  The  diffusion  rates  for  the  afternoon  model  show  a 
similar  behavior  to  the  prenoon  model  (Figure  5).  In  this 


model  <jJ,„  ~  0.45fte,  =  0.1 5Qe,  oj/c  =  0.05fte,  and  u>MC  = 

0.95fte.  The  peak  diffusion  rates  are  smaller  due  to  the 
lower  wave  amplitude  50  pT  and  are  higher  at  larger  pitch 
angles  than  the  prenoon  model  due  to  wave  power  extend¬ 
ing  to  lower  latitudes.  Clearly,  while  waves  in  this  sector 
contribute  to  loss  of  low-energy  electrons  and  acceleration, 
they  are  not  as  effective  as  the  prenoon  sector. 

5.  Global  Diffusion  Rates 

[22]  Energetic  electrons  drift  around  the  Earth  under  the 
influence  of  the  gradient  and  curvature  drift  and  the  con¬ 
vection  electric  field.  However,  for  high  energies  the  gra¬ 
dient  and  curvature  drift  are  the  most  important  factors.  The 
drift  periods  of  0.1  to  1  MeV  electrons  at  L  =  4.5  in  a  dipole 
magnetic  field  lie  between  approximately  100  and  16  min, 
respectively  Thus  during  a  storm  recovery  phase  lasting  a 
few  days,  electrons  complete  many  drift  orbits.  During  each 
drift  orbit  the  particles  will  encounter  waves  in  different 
MLT  sectors  and  be  scattered  in  pitch  angle  and  energy.  To 
take  this  into  account,  we  have  combined  the  results  of  the 
three  models  discussed  above  to  obtain  the  drift  average 
diffusion  rates.  We  have  multiplied  the  diffusion  rates  for 
each  model  by  the  fraction  of  MLT  that  the  model  repre¬ 
sents,  7/24  for  night,  6/24  for  prenoon,  and  3/24  for 
afternoon.  The  results  are  shown  in  Figure  6.  The  global 
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Figure  5.  Bounce-averaged  pitch  angle  (Daa)  and  energy 
(DEe)  diffusion  rates  for  whistler  mode  chorus  waves  as  a 
function  of  equatorial  pitch  angle  aeq  for  energies  of  30, 
100,  300,  and  1000  keV,  for  the  afternoon  model. 
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Figure  6.  Bounce-averaged  pitch  angle  and  energy  diffusion  rates  for  whistler  mode  chorus  waves  for 
the  night  (solid),  prenoon  (dotted),  and  afternoon  (dashed)  models,  weighted  according  to  the  MLT 
occurrence  of  the  waves.  The  energies  arc,  top  to  bottom,  1000,  300,  100,  and  30  keV. 


7  of  10 


A03225 


HORNE  ET  AL.:  TIMESCALE  FOR  ACCELERATION 


A03225 


diffusion  rate  due  to  waves  in  all  MLT  sectors  is  therefore 
the  sum  of  the  diffusion  rates  shown  in  each  panel  in 
Figure  6.  Thus  it  is  clear  that  energy  diffusion  can  extend 
to  large  pitch  angles  when  the  drift  of  the  particles  is  taken 
into  account. 

[23]  When  the  MLT  occurrence  of  the  waves  is  taken  into 
account,  only  waves  in  the  night  and  prcnoon  MLT  sectors 
contribute  significantly  to  particle  diffusion.  The  diffusion 
rates  generally  decrease  with  increasing  energy,  but  the  key 
factor  to  emerge  is  that  there  is  a  competition  between 
electron  acceleration  and  election  loss  that  is  energy  depen¬ 
dent.  By  assuming  a  separation  of  variables  and  ignoring 
mixed  pitch  angle-energy  diffusion,  the  timescale  for  elec¬ 
tron  loss  can  be  estimated  from  the  inverse  of  the  pitch 
angle  diffusion  rate  at  the  edge  of  the  loss  cone.  At  30  keV, 
the  loss  timescale  is  approximately  2  hours.  However,  at 
300  keV  the  loss  and  acceleration  timescales  are  comparable 
at  about  14  hours,  but  at  1  MeV  the  acceleration  timescale 
becomes  much  shorter  (3.9  days)  than  the  loss  timescale 
(23  days).  The  change  in  the  timescales  with  energy  is 
consistent  with  the  concept  of  energy  transfer,  that  as 
particles  are  diffused  into  the  loss  cone  at  low  energies  they 
give  energy  to  the  waves  which  accelerates  particles  at  high 
energies.  Since  the  timescale  for  loss  at  30  keV  is  a  few 
hours  and  since  ~30  keV  electrons  are  most  likely  to  be 
responsible  for  the  generation  of  chorus,  a  relatively  con¬ 
tinuous  source  of  ~30  keV  electrons  is  required  for  signif¬ 
icant  acceleration  up  to  a  few  MeV.  This  suggests  that 
prolonged  substorm  activity  is  a  requirement  for  efficient 
acceleration. 


Energy  (eV) 

Figure  7.  Time  evolution  of  the  increase  in  the  electron 
flux  due  to  whistler  mode  chorus  waves,  weighted  for 
latitude  and  the  occurrence  of  chorus  in  MLT.  The  flux  is 
shown  at  intervals  of  6  hours. 


F(E,  (xeq)  in  (9)  is  a  function  of  energy,  not  momentum,  and 
is  related  to  the  flux  J(Ey  ae?)  by 


‘V.O-tg/kO-JgSfrl'.*.).  (ID 


6.  Evolution  of  the  Particle  Flux 

[24]  The  timescale  estimates  calculated  above  are  for  the 
evolution  of  the  electron  distribution  function  and  not  the 
particle  flux  that  is  usually  measured  by  satellites.  To 
compare  with  satellite  observations,  we  first  calculate  the 
evolution  of  the  distribution  function  from  the  diffusion 
equation  and  then  convert  this  into  an  evolution  of  the 
electron  flux.  For  this  initial  assessment  we  neglect 
the  mixed  pitch  angle-energy  diffusion  coefficients  and 
use  the  rate  of  pitch  angle  diffusion  near  the  loss  cone  to 
calculate  the  timescale  for  losses  tl  to  the  atmosphere. 
Since  the  magnitude  of  the  mixed  diffusion  coefficients 
usually  lies  between  pure  pitch  angle  and  pure  energy 
diffusion,  the  results  are  most  likely  to  be  an  underestimate 
of  the  increase  in  flux.  Assuming  pitch  angle  isotropy,  the 
evolution  of  the  equatorial  energy  distribution  function  F(E , 
&eq)  can  be  obtained  by  bounce  averaging  (2)  and  writing 
(Dpp)  in  terms  of  the  bounce-averaged  energy  diffusion 
coefficient  (DE£)  and  is  given  by 


d_ 

OE 


A{E){DEE)ti{w)} 


F_ 

T  L  ’ 


(9) 


where 


A  =  (E  +  Eq)(E  +  2£o)^, 


(10) 


where  we  have  used /=  Jip2.  For  initial  conditions  we  use 
the  electron  flux  perpendicular  to  the  ambient  magnetic  field 
at  L-  4.5  measured  by  the  CRJRJES  satellite  on  13  October 
1990  (orbit  187)  during  the  recovery  phase  of  a  magnetic 
storm  [Meredith  et  al.,  2002a,  2002b].  The  flux,  measured 
at  18  logarithmically  spaced  intervals  between  0.28  and 
1.58  MeV,  was  interpolated  onto  a  high-resolution  regular 
grid  of  500  points  between  0.28  and  11  MeV  using  linear 
interpolation.  The  normalized  energy  diffusion  rates  ((O^)/ 
E2  =  5  x  10~5,  1 .5  x  1 0-5,  and  3  x  10~6  s_1)  at  energies  of 
100,  300,  and  1000  keV,  respectively,  were  also  interpolated 
onto  the  same  high-resolution  grid  together  with  the  loss 
timescales  obtained  from  the  pitch  angle  diffusion  rates  (rL- 
p2/{DaQ)  =  1.6  x  104,  5  x  104,  and  2.0  x  106  s).  An 
explicit  finite  difference  method,  accurate  to  second  order, 
was  used  to  solve  (9)  subject  to  the  condition  that  the  flux  at 
the  low-  (0.28  MeV)  and  high-energy  boundaries  (1 1  MeV) 
were  kept  constant. 

[25]  The  time  evolution  of  the  electron  flux  is  shown  in 
Figure  7.  In  the  absence  of  any  other  transport  processes, 
the  flux  >0.8  MeV  increases  by  more  than  an  order  of 
magnitude  and  approaches  a  steady  state  after  about 
24  hours  where  losses  balance  acceleration.  These  time- 
scales  are  comparable  to  the  observed  timescale  for  flux 
increase  in  the  radiation  belts  during  the  recovery  phase  of 
magnetic  storms  [Baker  et  al .,  1986,  1994;  Kim  and  Chan , 
1997;  Meredith  et  al ,  2002a,  2002b]. 


E  is  the  kinetic  energy,  E0  ~  m0c2  is  the  rest  mass  energy, 
and  where  the  last  term  represents  losses  to  the  atmosphere 
due  to  wave-particle  interactions.  The  distribution  function 


7.  Discussion 

[26]  Our  results  show  that  there  is  a  competition  between 
electron  acceleration  and  loss  for  diffusion  by  whistler  mode 
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chorus  waves.  Although  we  have  estimated  the  losses  from 
the  pitch  angle  diffusion  coefficients  at  the  edge  of  the  loss 
cone,  a  full  treatment  should  take  into  account  the  shape  of 
the  distribution  function  in  pitch  angle  and  energy.  The  pitch 
angle  diffusion  coefficients  are  relatively  flat  near  the  loss 
cone  and  are  well  below  the  strong  diffusion  limit.  There¬ 
fore  we  expect  that  the  pitch  angle  diffusion  rates  near  the 
edge  of  the  loss  cone  should  give  reasonable  estimates  of 
the  loss  rates. 

[27]  Although  losses  due  to  whistler  mode  chorus  are 
included,  losses  due  to  other  types  of  waves,  such  as 
plasmaspheric  hiss  and  electromagnetic  ion  cyclotron 
(EMIC)  waves,  have  been  omitted.  Plasmaspheric  hiss  is 
mainly  confined  to  within  the  high-density  plasmapause, 
which  is  compressed  during  large  geomagnetic  disturbances 
[Meredith  et  ai,  2004].  Since  wave  acceleration  is  most 
efficient  in  low-density  regions,  wave  acceleration  should 
take  place  outside  the  region  where  hiss  dominates.  Hiss 
should  contribute  to  radiation  belt  losses  during  the  later 
stages  of  the  recovery  phase  as  the  plasmapause  expands. 
Electron  losses  due  to  EMIC  waves  can  be  significant 
[Summers  and  Thorne ,  2003;  Albert ,  2003],  but  they  are 
generally  confined  to  electron  energies  greater  than  ~1  MeV 
[Meredith  et  ai9  2003b].  In  addition,  EMIC  waves  are 
usually  confined  in  MLT  to  regions  close  to  the  plasma¬ 
pause  boundary  and  therefore  electrons  would  only  be 
scattered  during  a  fraction  (<10%)  of  their  drift  orbit. 
Therefore  we  do  not  expect  these  omissions  to  impact 
significantly  on  the  main  result  presented  here. 

[28]  The  occurrence  of  whistler  mode  chorus  in  latitude 
and  MLT  may  have  a  significant  impact  on  the  acceleration 
of  radiation  belt  electrons.  Waves  in  the  prenoon  sector 
dominate  energy  diffusion  above  ~  1 00  keV.  Since  these 
waves  occur  at  latitudes  above  (and  below)  the  magnetic 
equator,  the  diffusion  caused  by  these  waves  peaks  at 
equatorial  pitch  angles  between  30°  and  50°.  If  pitch  angle 
diffusion  is  relatively  weak,  then  this  could  result  in 
butterfly  distributions.  The  results  presented  here  are  for 
conditions  averaged  over  the  lifetime  of  the  CRRES  orbit, 
but  one  may  envisage  particular  storms  where,  for  example, 
waves  at  high  latitudes  in  the  prenoon  sector  dominate 
waves  in  the  night  sector.  Under  these  conditions,  pitch 
angle  diffusion  would  be  very  small  for  pitch  angles  greater 
than  about  60°  and  butterfly  distributions  should  result. 
Flat-topped  pitch  angle  distributions  have  been  observed 
during  one  particular  magnetic  storm  during  the  recovery 
phase  when  the  electron  flux  was  observed  to  increase 
[Horne  et  a /.,  2003a].  These  distributions  were  observed 
near  L  *=  4.5,  where  magnetopause  shadowing  can  be  ruled 
out  and  where  drift  shell  splitting  should  not  be  significant. 
We  suggest  that  in  addition  to  drift  shell  splitting,  flat- 
topped  distributions  and  butterfly  distributions  could  be 
produced  as  a  result  of  wave  acceleration  and  could  be  a 
signature  of  wave  acceleration.  A  careful  analysis  of  particle 
pitch  angle  distributions  in  MLT  is  required  to  test  this 
suggestion. 

8.  Conclusions 

[29]  We  have  used  wave  and  particle  data  from  the  CRRES 
satellite  to  determine  the  timescale  for  electron  acceleration 
in  the  heart  of  the  outer  radiation  belt  near  L  -  4.5  due  to 


whistler  mode  chorus  waves.  During  magnetically  active 
conditions,  we  split  the  wave  occurrence  into  three  MLT 
sectors  and  determined  the  contribution  to  wave  acceleration 
and  loss  from  each  sector  by  calculating  the  bounce  average 
pitch  angle  and  energy  diffusion  rates.  The  results  were 
weighted  by  the  occurrence  of  the  waves  in  MLT  to  obtain 
the  timescale  for  wave  acceleration.  The  main  conclusions 
are  as  follows: 

[30]  1.  Chorus  wave  amplitudes  in  the  night  MLT  sector 
peak  near  the  magnetic  equator  and  contribute  to  significant 
pitch  angle  and  energy  diffusion  at  pitch  angles  between  0° 
and  85°.  However,  in  the  prenoon  and  afternoon  MLT 
sectors  chorus  wave  amplitudes  peak  at  latitudes  above 
and  below  the  magnetic  equator.  These  waves  are  only 
effective  in  causing  pitch  angle  and  energy  diffusion  at 
equatorial  pitch  angles  less  than  60°. 

[31]  2.  At  energies  less  than  ~300  keV,  there  is  a 
competition  between  electron  acceleration  and  loss  to  the 
atmosphere.  Above  ~300  keV,  electrons  are  accelerated 
faster  than  they  are  lost. 

[32]  3.  After  bounce-averaging  and  weighting  the  dif¬ 
fusion  coefficients  by  the  occurrence  of  whistler  mode 
chorus  in  MLT,  we  find  that  the  timescale  for  an  order 
of  magnitude  increase  in  the  electron  flux  at  l  MeV  is 
about  24  hours  at  L  =  4.5.  This  compares  favorably  with 
observations. 

[33]  4.  We  suggest  that  electron  acceleration  by  whistler 
mode  chorus  waves  on  the  dayside  of  the  magnetosphere 
could  produce  butterfly  pitch  angle  distributions.  Since  drift 
shell  splitting  becomes  small  at  low  L,  the  best  place  to 
determine  whether  wave  acceleration  could  be  responsible 
for  butterfly  distributions  is  at  low  L  during  large  geomag¬ 
netic  storms. 

[34]  We  conclude  that  wave  acceleration  by  whistler 
mode  chorus  is  a  viable  mechanism  to  explain  electron  flux 
increases  during  magnetically  disturbed  times  in  the  heart  of 
the  radiation  belts  near  L  =  4.5. 
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